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Ground State Phase Diagram of the One Dimensional S — 1/2 XXZ 
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The one dimensional S — 1/2 XXZ model with dimerization (1 — j) and quadrumerization 
8 is studied by the numerical exact diagonalization of finite size systems. Using the conformal 
field theory and the level spectroscopy method, we calculate the ground state phase diagram 
with XY-like anisotropy A (0 < A < 1). The ground states of this model contain the Haldane 
state, S — 1 dimer state, S = 1 large- D state and 5 = 1/2 dimer state as limiting cases. The A 
and ^-dependence of the critical exponent v of the energy gap is determined from the conformal 
dimensions of excited states. 
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§1. Introduction 

Recently, the one dimensional S = 1/2 Heisenberg 
chains with a variety of spatial structures have attracted 
a great deal of attention. Although the uniform S = 1 12 
XXZ model can be solved by the Bethe's hypothesise! 
and is known to have a gapless ground state for the XY- 
like anisotropy, the various spatial structures such as the 
spin-Peierls dimerization and ladder structure can drive 
this ground state into thCj-spin gap states. 

In our previous work,tf the phase diagram of the 
S = 1/2 isotropic Heisenberg model with dimerization 
and quadrumerization was studied by the exact diago- 
nalization and phenomenolozical renormalization group 
method. The XY case is also studied analytically. How- 
ever, the phase diagrams in these two limiting cases are 
quite different. Therefore it is worthwhile to investigate 
the crossover between these two limiting cases. 

This paper is organized as follows. In the next sec- 
tion, the model Hamiltonian is defined. The numerical 
results are presented in §3. Analyzing the exact diag- 
onalization data by the conformal field theory and the 
level spectroscopy method, we calculate the exponent v 
characterizing the opening of the energy gap. The last 
section is devoted to summary and discussion. 

§2. Model Hamiltonian 

The Hamiltonian of the one dimensional dimerized and 
quadrumerized 5=1/2 XXZ chain is given by 
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where 1 — j(— oo < j < oo) and 5(— 1 < S < 1) rep- 
resent the degree of dimerization and quadrumerization, 
respectively. The anisotropy parameter is denoted by A. 
In our previous work,EP we calculated the phase diagrams 
for the cases A = 1 (isotropic case) and (XY case). In 
the present work, we concentrate on the XY-like regime 
< A < 1. 

§3. Numerical Results 

3. 1 Phase diagram for < A < 1 

In order to determine the phase boundary with high 
accuracy, we use the .twisted boundary method of Ki- 
tazawa and Nomura.Q'El' The Hamiltonian is numerically 
diagonalized to calculate the two low lying energy lev- 
els with the twisted boundary condition S^ N+1 = —Sf, 
SIn+1 = Si S% N+1 = Sf for AN = 12, 16, 20 and 24. 

From the result of the isotropic and XY cases, El 1 we 
expect two different kinds of ground states. For small 
values of j and 6, the ground state is the Haldane-like 
state with valence bond solid (VBS) structure. Under 
the twisted boundary condition, the eigenvalues of the 
space inversion P and the spin reversal T are all equal 
to — loW As j and/or 8 increases, the phase transition 
takes place into the dimer-like state for which P = 1 
and T = 1. We make use of P and T eigenvalues to 
distinguish the two phases. For fixed 8, the energies of 
the two states vary with j. For small | j |, the energy of 
the Haldane state is lower than that of the dimer state. 
As | j | increases, the latter becomes lower than the 
former. The two levels cross at one point which gives the 
finite size transition point j = j c (N). Fig. 1 shows the 
^-dependence of the two lowest levels for N = 24, A = 
0.5, 8 = 0.4 and j > 0. We extrapolate the critical point 
as jc(N) = j c (oo) + aN- 2 + bN- i . 

In Fig. 2, the ground state phase boundaries are rep- 
resented by the solid line, A, □ and o for A = 0,0.5,0.7 
and 1, respectively. For < A < 1, the phase boundaries 
are closed at the negative finite value of j = j c (8) . This 
is in contrast to the isotropic case in which the Haldane 
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phase remains stable even in the limit j — ► — oo as far as 
| 8 |< 5 C — 0.25. Here the critical value 8 C corresponds to 
the critical dimerization of the Jlaldane-dimer transition 
of the 5 = 1 Heisenberg chain.BQ' This behavior can be 
understood in the similar way as the XY case© and the 
ground state for j < j c (8) corresponds to the large- D 
phase of the 5=1 Heisenberg chain for | 8 |< 8 C . It 
should be noted that there is no critical point between 
the large- D phase and the dimer phase_as in the case of 
the 5=1 dimerized Heisenberg chain.® 

For j ^ and 8 — ► 1, the spins connected by the 1 + 8- 
bonds form strong singlet pairs. The effective coupling 
j e jf and effective anisotropy A e g between the spins Su+i 
and can be calculated by the perturbation theory 

as, 
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Therefore the effective Hamiltonian is given by 
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We have also diagonalized this effective model to check 
the phase boundary for small j. The phase boundary 
obtained from the effective model is shown by x in Fig. 
2 for the case A = 0.5. 

3. 2 Critical exponent v of the energy gap 

At the critical points, we calculate the energy with 
fixed wave number k by Lanczos exact diagonalization 
method under the periodic boundary condition for the 

4N 

magnetization M z = and 1 where M z = 5f . The 

2=1 

system sizes are AN — 12, 16, 20 and 24. The ground 
state has M z = and k = 0. 

It is known that the finite size correction to the ground 
state energy is related with the .ceiijbpLl charge c and the 
spin wave velocity v s as follows p t3 
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where E g (N) is the ground state energy, Ek 1 (N), the 
energy of the excited state with k\ = ^ and M z = 0, c, 
the central charge and £oo, the ground state energy per 
unit cell in the thermody nam ic limit. After appropriate 
extrapolation to N — > oo JEB we have checked that the 
central charge c is close to unity on the phase boundary. 
Therefore, the present model can be described by the 
Gaussian model on the critical line, 
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where <f> is the boson field defined in the interval < 
(j) < \p2m and n is the momentum density conjugate to 
4> which satisfies [4>(x),Il(x')] = i8(x— x'). The deviation 
from the critical point is described by the term cos y2<f> 
and the low energy properties of our model can be de- 
scribed by the one-dimensional quantum sine-Gordon 
theory near the critical line, 
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The scaling dimensions x n of the operators are related 
with the energy eigenvalues E n (N, M z ) of the corre- 
sponding excited states as x n = limjv^oo x n (N) where, 

N 

Xn{N) = [E n (N,M z ) - Eg(N)]. (3.8) 
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We can identify the correspondence between the op- 
erators in boson representation and the eigenstates_af 
spin chains by comparing their symmetry properties .trE3 
Let us denote the scaling dimensions of the operators 
v / 2cosv / 20 and v^sinv^^ by x\ and X2, respectively. 
Both of them should be equal to K/2 in the thermody- 
namic limit as determined from their correlation func- 
tions. Numerically, these exponents are calculated using 
Eq.Q for M z = 0, P = 1, k = 0[xi] and M z = 0, 
P = —1, k = 0[x2]. Thus the value of K can be deter- 
mined from these values. Actually, to reduce the finite 
size correction to 0(1/A^ 2 ), it is more convenient to use 
the combination, 

K(N) = Xl {N) + x 2 (N), (3.9) 

as proposed by Kitazawa and Nomura.!! 

For j > 0, the phase boundary approaches that of 
the isotropic chain as A tends to unity and the logarith- 
mic corrections appear due to the SU(2) symmetry of 
the problem. For j < 0, this happens only for 8 >■ A,. 
According to the renormalization group calculationptiP 
these logarithmic corrections are given by, 
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where x^ is the scaling dimension of the operator 
exp(±i#), yoil) is proportional to 1/logiV at the SU(2) 
symmetric point and the deviation from the SU(2) sym- 
metric point is denoted by t. The scaling dimension 24 
is calculated by Eq.(O) for M z = 1, P = 1, k = 0. In 



this case, we use the combination 
1 



K(N) 



K(N) 



= 21 (AO +22 (AO + 2x4 (AO, (3-13) 



to reduce the influence of the logarithmic corrections for 
A - 1. 
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In both cases, we assume the formula 



A-(.V) + ^ + 



(3.14) 



to extrapolate K(N) to N — > oo. 
v of the energy gap is given by 

1 
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The critical exponent 
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For A > 0.9, we use eq. (3.13) as far as j > or 
j < and 5 > <5 C . In other cases, eq. ( fTs| ) is used. The 
extrapolation procedure of 7^ to N — > oo is shown in Fig. 
3 for A = 0.9, S = 0.4 and j = 0.9354. Figure 4 shows 
the A and S dependence of the critical exponent v. The 
error bars are estimated from the difference between the 
values extrapolated from 12 < AN < 20 and those from 
12 < AN < 24. The values of K estimated from other 
combinations coincide with the above estimation within 
the error bars. As shown in Fig. 4 the critical exponent v 
is almost independent of 5 and decreases with increasing 
A for j > 0. At A = 1 and 0, the critical exponent 
v reduces to 2/3 and 1, respectively. For S = and 
j = 1, the parameter K can be calculated by the Bethe's 
hypothesisESEiP as 

(3.16) 



K = 



arccosA 

The numerical results are in agreement with the exact 
results as shown in Fig. 4. For j < 0, the A-dependence 
is rather complicated. For | 6 |> S c , the phase bound- 
ary approaches that of the isotropic model as A tends 
to unity. At the same time, the critical exponent v ap- 
proaches 2/3 which is the value for the SU(2) symmetric 
critical point. On the other hand, for | S |< <5 C , the 
phase boundary tends to j — > — oo as A — > 1. As dis- 
cussed in the preceding section, the phase transition in 
this limit corresponds to the Haldane-large-D transition 
in the 5 = 1 dimerized XXZ chain. Therefore, the criti- 
cal exponent v approaches the nonuniversal value which 
varies with S. Although the precise estimation of v for 
the S = 1 chain is not available in literature, Glaus and 
SchneiderEZP estimated as v ~ 1.5 at S — 0. Our results 
are consistent with their estimation. 

Our value of v for S = is also consistent with the 
critical exnonents obtained by Yamanaka, Hatsugai and 
KohmotolHr and Yamanaka and Kohmoto.t£l ) It should 



it the sign of exponent a is inverted in Fig. 
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§4. Summary and Discussion 

The ground state phase diagram and the critical ex- 
ponent v of the dimerized and quadrumerized spin-1/2 
XXZ chain is calculated by the numerical diagonaliza- 
tion method for < A < 1. The critical points are de- 
termined by the method of twisted boundary condition. 
For 5 ~ 1 and j ~ 0, the numerical results are compared 
with those for the effective model obtained perturbation- 
ally. The gap exponent v is calculated by analyzing the 
numerical diagonalization data using the conforma! field 
theory. 

For j > 0, the phase boundaries are insensitive to A, 



while for j < 0, they are closed at finite negative values 
of j = j c (S) and the large- D phase appears in the ground 
state for j < j c (S). 

The A-dependence of the exponent v is qualitatively 
different according as j > or j < 0. For j > 0, the crit- 
ical exponent v is insensitive to S and always reduces to 
the SU(2) symmetric value 2/3 as A — > 1. For j < and 
S < S c , v approaches the nonuniversal value which corre- 
sponds to the value for the large- D-Haldane transition in 
the S = 1 dimerized XXZ model as A approaches unity. 
On the other hand, it approaches the SU(2) symmetric 
value 2/3 at A = 1 for j < and 6 > S c . 

In this paper, we concentrated on the parameter 
regime — 1 < S < 1 and < A < 1 to clarify the feature 
of the crossover between the XY case and the isotropic 
case. It is expected, however, that this model has more 
variety of phases outside this parameter regime— -Espe- 
cially, from the consideration of the case 5 = 0,E£lE3' we 
may expect the presence of the XY and ferromagnetic 
phase for A < and the Ising phase for A > 1. It is 
interesting to find how these phases are modified in the 
presence of the quadrumerization. The investigation of 
these problems is left for future studies. 

We thank A. Kitazawa, K. Okamoto, S. Yamaguchi 
and M. Yamanaka for valuable comments and discus- 
sion. We are also grateful to H. Nishimori for the pro- 
gram package TITPACK version 2 for the diagonaliza- 
tion of spin-1/2 systems. The numerical calculation is 
performed using the HITAC S820 and SR2201 at the In- 
formation Processing Center of Saitama University and 
the FACOM vpp500 at the Supercomputer Center of In- 
stitute for Solid State Physics, University of Tokyo. 
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Fig. 1. The j-dependence of the two lowest energy eigenvalues 
with twist boundary condition. The energies of the Haldane 
state and dimer state are represented by o and •, respectively, 
for N = 24, A = 0.5, 6 = 0.4 and j > 0. 



Fig. 2. The phase diagram of the isotropic model (o), A = 0.7 
(□), A = 0.5 (A) and the XY model (solid line). The result for 
the effective model (3.3) is shown by X for A = 0.5. The broken 
lines are the guide for the eye. 



Fig. 3. The extrapolation procedure of finite size K for A = 0.9, 
5 = 0.4 and j > 0. 



Fig. 4. The A dependence of the numerically obtained critical 
exponent v. The open and filled symbols represent the cases 
j > and j < 0, respectively. For j > 0, the 5-dependence is 
almost unvisible. 
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